Entanglement entropy in a boundary impurity model 
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Boundary impurities are known to dramatically alter certain bulk properties of 1 + 1 dimensional 
strongly correlated systems. The entanglement entropy of a zero temperature Luttinger liquid 
bisected by a single impurity is computed using a novel finite size scaling/bosonization scheme. For 
a Luttinger liquid of length 2L and UV cut off e, the boundary impurity correction (<5Si mp ) to the 
bulk logarithmic entanglement entropy (Sent oc InL/e) scales as SSimp ~ y r \nL/e, where y r is the 
renormalized backscattering coupling constant. In this way, bulk entanglement entropy within a 
region is related to scattering through the region's boundary. In the repulsive case (g < 1), SSimp 
diverges (negatively) suggesting that the bulk entropy vanishes. Our results are consistent with 
the recent conjecture that entanglement entropy decreases irreversibly along renormalization group 
flow. 
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Quantum field theories describe coupled quantum os- 
cillators and, therefore, even the ground state of nonin- 
teracting fields may exhibit quantum correlations over a 
long range. Consider the ground state of the free electro- 
magnetic field. An electric field amplitude E at a point x 
is quantum correlated with an amplitude E' at a point x' 
in that there are contributions to the ground state wave- 
function containing such products as . . . \E) x \E) x i . . .. 
Since the ground state contains superpositions of such 
states, it is said to be entangled — even though it is a free 
field. If, in a d-dimensional spacetime, a spatial box of 
volume is formed, it follows that the degrees of free- 
dom which reside exclusively in the box will appear to be 
in a mixed state. The degree of mixing may be character- 
ized by the von Neumann entropy, S = — trplnp, where 
the reduced density matrix p has been formed by tracing 
over the degrees of freedom exterior to the volume L 1 . 

Geometric or entanglement entropy formed in this 
fashion was introduced in the context of black hole 
quantum mechanics and Hawking-Bekenstein entropy , 
where it was found that entanglement entropy is not 
an extensive quantity but, rather, scales as the area of 
the bounding surface, S cx L d ~ 2 2J. This highly sug- 
gestive result is believed to bear some relation to holo- 
graphic principle proposals In condensed matter 
physics, the role of entanglement entropy in understand- 
ing quantum critical phenomena 4] and quantum phase 
transitions has recently been emphasized. 1 + 1 di- 
mensional conformal field theories — which describe crit- 
ical spin chains, Luttinger liquids and other massless 
theories — have pointlike bounding surfaces; remarkably, 
the entanglement entropy was shown to depend uni- 
versally upon the central charge of the theory and to 
diverge logarithmically with the length of the subsys- 
tem ja, H HI Specifically, the entropy is given by 
S — | InL/e where c is the central charge. 

Studies of entanglement entropy have, so far, been re- 
stricted to homogeneous models. [For a recent excep- 
tion, see E3.] However, consider the artificial introduc- 



tion of an inhomogeneity at one boundary point. Quan- 
tum entanglements across the boundary must be associ- 
ated with scattering through the boundary and, there- 
fore, a boundary impurity that alters the bulk conduc- 
tance is expected to alter the bulk entanglement entropy. 
Specifically, backscattering from an impurity should be 
associated with negative corrections to the bulk entan- 
glement entropy — if a sensible perturbative approach to 
this highly non-local quantity can be found. In this 
manuscript, we develop a scheme by which these correc- 
tions to entropy may be computed for a general bound- 
ary impurity model and demonstrate a connection be- 
tween the scaling of backscattering and that of entangle- 
ment entropy for the specific model of an impurity in a 
Luttinger liquid (LL). Significantly, we find a logarith- 
mic divergence in the entropy computed in this fashion 
which has its origin in an infrared divergence in the effec- 
tive boundary theory. Finally, our perturbation theory 
directly addresses and supports the conjecture that 
entanglement entropy decreases irreversibly along renor- 
malization group flows. 

We adopt the model of Kane and Fisher (KF) [irT |. 
in which a single impurity in a Luttinger liquid was 
shown to have a dramatic effect — either effectively de- 
coupling the two sides, or effectively vanishing, depend- 
ing on whether the LL is repulsive (g < 1) or attractive 
(g > 1). The behavior of £ e nt suggests the following 
puzzle: If one considers (following |£|) a "mixed" real 
space/momentum space basis where at wavenumber k 
there are a series of (kL) spatial modes of width k^ 1 , 
entanglement entropy arises from those spatial modes 
which "straddle" the boundar y; c ounting them naively 
gives In (A: max /A: m i n ) ~ InL/e The UV divergence 

of the entropy is a consequence of the unlimited number 
modes close the boundary. If an impurity is introduced 
into a repulsive LL, the effective strength grows with 
inverse wavenumber and therefore the contributions to 
entanglement become weak at the longest wavelengths. 
Using the KF scaling of impurity strength, one might es- 
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timate the correction to entanglement entropy, 8Si mp , by 
excluding contributions to the sum above from wavenum- 
bers less than V^L/e) 1 ^ 9 , the approximate energy scale 
of the impurity in a finite lattice (1/T ^S> L). This esti- 
mate gives JS'imp ~ (g — 1) In (L/e) — const which demon- 
strates the expected reduction. However, this estimate 
suggests that S cn t may still diverge as L/e — > oo, contrary 
to the expectation that the impurity effectively "discon- 
nects" the two sides for g < 1. 

To proceed, Sent must be directly calculated or ex- 
pressed as operators with known LL scaling dimensions. 
However, S en t is an unusual, non-local quantity and de- 
spite its superficial resemblance to a log-scaling marginal 
operator, it has no simple operator representation within 
the LL. An elegant geometric approach to this problem 
was pioneered by Callan and Wilczek Q and Cardy ^3 > 
who introduced a euclidean path integral representation 
for the geometric entropy of a 1 + 1 dimensional field the- 
ory computed on the half-infinite spatial domain. First 
expressing the ground state wavefunctional as a euclidean 
path integral 

V[(f>(x)] = (0|O) cx J D(f>e- S ^\ (1) 

the ground state density matrix functional may be writ- 
ten as a double path integral where one identifies the 
imaginary time domain to be r € (—oo, 0) for one path 
and r £ (oo,0) for the other. Now the reduced density 
matrix functional for the fields 4>± on halffine (x > 0), 
(4>+\p\4>-) , may be expressed as a path integral over the 
entire 1 + 1 dimensional euclidean space (x, t 6 (— oo, oo)) 
with <f>± as boundary conditions on cither side of a cut 
along the polar angle 9 — (figure 1). If this cut is en- 
larged to a finite "deficit angle," a, the path integral is 
now evaluated on a cone. 

Z a = (0+M0_) = / Dcpe- 3 ^ (2) 

J Ft 2 / {6=0} 

Using the replica trick, it was shown that the entropy 
Sent = — trplnp may be written as a variation of the 
partition function, Z a , with respect to a 0. Specif- 
ically, Sent = ( — 27r^j + 1) lnZ a \ a —Q. Remarkably, this 
expression has the same form as thermodynamic entropy, 
a = (— fl-Mj + 1) hi Z, where /3 is identified with the com- 
plete polar angle 2tt. 

The action appearing in @ is taken to be the LL ac- 
tion, So, plus a weak impurity potential scattering term, 
Simp. To express the boundary conditions imposed in 
equation the action is written in polar coordinates: 

S a = - [ dff f rdrhd r <t>f + \(d e ^) 2 ] (3) 

9 JO J a 1 r 

+ y / dr cos V^TT<p(r, 9) (4) 

e=±f 




FIG. 1: Conical geometry corresponding to path integral 
Deficit angle, a, shown here is negative. Dotted line is the 
(spacelike) impurity worldline. 

Without the impurity, the computation of S cn t, at zero 
physical temperature, amounts to the computation of a 
conventional thermodynamic entropy from the free en- 
ergy but on a conical manifold with the (imaginary) time- 
like angle variable having period [3 = 27r + a. The cor- 
rections to the entropy from the boundary, SS- lmp , then 
involve the finite temperature (f3) corrections of the im- 
purity free energy, SF, to the bulk free energy, F; specifi- 
cally, the variation of 5F with respect to a about [5 = 27r. 

SS imp = (-2n^- + l)(—(35F)\ a=0 (5) 



However, as depicted in figure 1, the impurity world- 
line has become spacelike. This state of affairs is, in some 
sense (see figure 2), analogous to the computation of fi- 
nite length corrections to the pressure in a Luttinger liq- 
uid due to a single impurity. The correction to the pres- 
sure is a variation of free energy with respect to length: 
Sp = —dSF/dL, where 5F is the lowest order correction 
to the bulk free energy coming from the impurity. In the 
finite size scaling regime in which L <C (3, SF, and there- 
fore, Sp, will depend algebraically on length, although 
they both vanish thermodynamically (i.e. SF ~ L 1 ^ 29 ; 
Sp ~ L~ 2g ). In the entanglement entropy problem, the 
finite temperature scaling regime is reached when the 
quantization length scale exceeds the inverse tempera- 
ture, (3 — 27r. It will be shown that this condition is 
In (L/e) ^> (3, which can always be met. 

To compute SF, the degrees of freedom in the bulk 
will be integrated out to produce an effective boundary 
sine-Gordon action for the "past" and "future" impurity 
fields at 9 = ±tt/2. There are two complicating charac- 
teristics: 1) The two impurity fields interact and 2) The 
impurity action lacks translation invariancc which inter- 
feres with the standard momentum shell rcnormalization 
group scheme. 

The free action is brought into a separable gaussian 
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FIG. 2: Analogous computation of finite L correction to the 
pressure in the presence of an impurity (dotted line). 



form with the "log-periodic" expansion: 

1 = ^ sin (k m In r 

t 



sin(fc ro ln-)e- <w "V(fe,w) (6) 



where we introduce Matsubara frequencies cu n = =|p = 
i+a/2-K * na * explicitly incorporate the deficit angle, 
wavenumber k m = ^S- and P = In The spatial modes 
are made to vanish at L and at a radial UV cut off, e. 
The action now becomes: 



S 



(7) 



Computation of 5 cn t proceeds exactly as the standard 
thermodynamics computation except the final replace- 
ment of the sum over k m by an integral now includes a 
logarithmic density of states, lni/e. Within the (Cardy- 
Callan-Wilzcek) conical geometry scheme, this is origin 
of the logarithmic entanglement entropy. 

Following the standard procedure for deriving the ef- 
fective impurity action, we introduce a pair of constraint 
fields, A±(fc), and impurity fields, 4>±(r), which corre- 
spond to the past and future branches of the impurity 
worldline. After the bulk fields (j)(rff) are integrated out, 
the action for the constraint fields is: 



Sc. = 





(8) 



k m ;a+± 



where d(k) = k m tanh . We have also made use of the 
fact that (f>*(ku>) = <j>(k — ui), etc. Finally, we integrate 
out the constraint fields and diagonalize the action using 
cf>± = -j^{4> ± </>)■ The final effective impurity action is 
then: 
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Sqi + Si 
- d(k m ) 



2y / dr cos 
Jo 



(9) 




1+A(fe m ) 






r) cos v 2n(j)(r) 



where A(fc m ) = l/cosh^ 22 -. This impurity model 
closely resembles the double barrier impurity model of 
[I7)j | although the impurity worldlines are spacelike. Sim- 
ilarly, the field <p, represents the number of particles "be- 
tween" the two barriers and has an energy gap of 0(1/ (3). 
There is no comparable "gate voltage" in this problem 
and thus (cf>) = 0. Only the mode <f> is active and we 
reach the final effective low energy action: 



S = -^<i(fc m )0 2 (fc m ) + y drt 

3 u ** e 



cos V47r0(r) (10) 



As expected, the finite time scale, (3, appears in the ef- 
fective spacelike model for the impurity. In an ordinary 
impurity model, renomalization group flow would termi- 
nate at the time scale and thus the algebraic behavior 
free energy could be directly determined. However, re- 
calling the "log-periodic" substitution, equation JBJ, this 
model no longer has translation invariance in r. For this 
reason, it is not obvious how to implement a momen- 
tum shell renormalization group calculation. Instead, we 
directly compute the corrections to the bulk free energy 
arising from the impurity potential in (|10|l to lowest order 
in y: 



(35F^ = ]y 2 J^ dnj dr 2 G(n,r 2 ) 



where 



G(n,r 2 ) 



— I i\fiTTtj>(ri) -i\/Tn4>(r2)\ 



(11) 



(12) 



G(ri,r 2 ) is most easily calculated in terms of logarithmic 
variables, p = lnr/e: 



G{pi,P2) = 



(sinh ^+^) )29( s inh *(Pi-P'> )-2fl 



(sinh^)s(sinh^p)s 



(13) 



From the form of G(p\p 2 ) 1 the free energy correction 
8F^ will only be in the finite (3 scaling range (and thus 
depend upon (3 and the deficit angle, a, as desired) if the 
characteristic logarithmic length scale, p 3> (3- Therefore, 
we arrive at the condition for finite (3 scaling: InL/e 3> [3. 

Returning to radial variables, the "log-periodic" be- 
havior is ultimately seen to leads to a large radial dis- 
placement Green's function that is algebraic in the ratio 
of radii, 



G(r ir2 ) = 



-29 



(14) 



The IR divergent part of the correction 5F l|ll|) comes 

g 

from r 2 ~ e and r\ ~ L where G ~ (r 2 /ri) 1+a / 2w and 
thus we arrive at the leading order divergence in the im- 
purity correction to the free energy: 



[35F 



(2) = 
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1 - (2^g/(3f 
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(15) 



Finally, the boundary impurity correction, SSi mp , to the 
bulk entanglement entropy is found from the variation of 
SF (equations © and l|TH)l. for g ^ 1): 



imp 



y 2 e 2 



L 



1-9 



In - - c(g) 



(16) 



where c{g) is a positive constant, independent of L/e. 
For 5 ->• 1, 



y 2 e 2 



(l+.9) : 



g In — 



(17) 



For L/e 1, the entropy correction is strictly nega- 
tive for all values of g. The appearance of the inverse 
temperature, (3, (or equivalently, the deficit angle a) in 
the exponent of the free energy 1)15(1 is responsible for the 
negative value of the entropy, (JS'imp, independent of g. 
In a free massless gas, the free energy increases to zero 
as the temperature goes to zero (or /3 — > oo) and thus 
the variation (0 yields a positive entropy. Here, the free 
energy decreases as f3 — > 27r~. 

For repulsive Luttinger liquids (g < 1), 5Sj mp diverges 
algebraically in L/e suggesting that the bulk entangle- 
ment entropy, in the presence of an impurity, is zero for 
sufficiently large L. Note, however, that for finite (bare) 
coupling y and finite L/e, the entanglement entropy is 
nonzero, suggesting that high energy modes close to the 
boundary remain entangled until a critical coupling or 
length scale is reached. For attractive Luttinger liquids 
(g > 1), SSi m p approaches a negative constant with alge- 
braic corrections in L/e that vanish as L/e — > oo; there- 
fore, an attractive LL remains maximally entangled. For 
the noninteracting case (g = 1), the correction diverges 
as In 2 L/e suggesting that, even for the noninteracting 
fermions, the boundary impurity suppresses the entan- 
glement entropy. 

Defining j/o = V e as the dimensionless coupling on the 
scale of the lattice and y r = yo (L/e) 1-9 as the renor- 
malized coupling on the scale, L (using the RG results 
of 0), SSi mp = -yoj/rln-jf- Our perturbative results 
may then be improved by the renormalization group and 
written in the following suggestive way: 



Sent remains logarithmic, but with a universal prefactor, 
c which is smaller than the central charge of the clean 
system, consistent with this proposal. In the zero physi- 
cal temperature KF impurity model, RG flow is accessed 
through L/e — > oo (which also drives the system further 
into the finite (3 scaling limit on which our results are 
based.) Thus, noting the sign change in the ^-dependent 
prefactor of ifTfljl . SSi mp strictly decreases for all g along 
the renormalization group trajectory set by yo and in- 
creasing L/e. 

Lastly, we find it significant that the logarithmic be- 
havior of the entropy (eqns. ifTfijl . ifTTjl and JIHJ) is now 
a symptom of an IR divergence in the boundary theory, 
rather than the logarithmic divergence of the wavenum- 
ber eigenvalues (fc" 1 oc In L/e) responsible for the bulk 
entanglement entropy within the Cardy-Callan-Wilczek 
construction. 

It should be noted that these results pertain to a LL 
bisected by an impurity; it is unknown whether this re- 
sult remains true for a region of length L imbedded in an 
infinite LL. (Recently, Calabrese and Cardy have given 
generalizations to the bulk entanglement entropy for un- 
equal partitions of 1 + 1 dimensional CFTs Q.) Further- 
more equation (|18f) does not reflect the contribution to 
the entropy from the open boundaries at ±L implicit in 
our construction [tj. 

The bosonization scheme presented here may likely be 
applied to Kondo type impurities as well as spin chains 
with comparable boundary impurities. Entanglement of 
an impurity spin with conduction channels plays a crucial 
role in more complicated boundary problems such as the 
two-channel Kondo effect. Furthermore, the predictions 
for the behavior of the entanglement entropy in the pres- 
ence of an impurity might be checked by density matrix 
renormalization group (DMRG) methods, since DMRG 
explicitly follows the flow of eigenvalues of the reduced 
density matrix. 
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Alex Zaharakis for Mathematica support, and the Kavli 
Institute for Theoretical Physics at UCSB where this 
work was completed; this research was supported in part 
by the National Science Foundation under Grant No. 
PHY99-07949. 
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for g < 1 (using c = 1 for a spinless LL). In |4J, en- 
tanglement loss was suggested as a fundamental fea- 
ture of the irreversibility of renormalization group flow. 
This intriguing proposal followed alon g th e lines of the 
"c-theorem" of conformal field theory |l3j |. Refael and 
Moore 01 have recently computed the entanglement en- 
tropy for strongly disordered spin chains and found that 



[1] N. D. Birrell and P. C. W. Davies, Quantum Fields in 
Curved Space, Cambridge University Press (1982). 

[2] M. Srednicki, Phys. Rev. Lett. 71 666 (1993) 
[hep-th/9 303048 . 

[3] R. Bousso, Rev. Mod. Phys. 74, 825 (2002) and refer- 
ences therein. 

[4] G. Vidal, J. I. Latorre, E. Rico and A. Kitaev, Phys. 
Rev. Lett. 90, 227902 (2003); J. I. Latorre, E. Rico, 
and G. Vidal, Quant. Inf. and Comp. 4, 048 (2004) 



5 



quant-ph/03040 98|; J. I. Latorre, C. A. Lutken, E. Rico, 

and G. Vidal, |quaiirp h/0404120]. 
[5] A. Osterloh, L. Amico, G. Falci and R. Fazio, Nature 416 

608 (2002); T. J. O sborne and M. A. N ielsen, Phys. Rev. 

A 66 32110 (2002) quant-ph/0202162 ; F. Verstraete, M. 

A. Martin-Del gado, and J. I Gira c, Phys. Rev. Lett. 92 

087201 (2004) quant-ph/0311087 . 
[6] C. G. Callan and F. Wilczek, Phys. Lett. 333B, 55 (1994) 

hep-th/9401072 . 
[7] C. Holzhe y, F. Tarsen, and F. Wilczek, Nucl. Phys. B424 

44 (1994) hep-th/9403108 . 
[8] T. M. Fiola, J. Preskill, A. Strominger, and S. P. Trivedi, 

Phys. Rev. D50 hep-th/9403137]. 



[9] P. Calabrese and J. Cardy hep-th/0405152 . 
[10] C. L. Kane and Mathew P. A. Fisher, Phys. Rev. B 46, 
15233 (1992). 

[11] At scale n, the wavenumber k — 2 n fc m m contributes 
O(kL) spatial modes, only one of which straddles the 
boundary and is entangled. However, the maximum num- 
ber of "binary" scales is n max ~ log fc max /fc m i n , where 
fcmax is the UV cut off e _1 . 

[12] J. Cardy and I. Peschel, Nucl. Phys. B300 377 (1988). 

[13] A. B. Zamolodchikov, JETP Lett. 43, 730 (1986). 

[14] G. Refael and J. E. Moore, cond-mat/0406737 . 



